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INTRODUCTION 


Essentially all of the important structural problems related to the 
design of high-temperature energy system components, e.g., aircraft and 
rocket engine components, and nuclear reactor system components, are non- 
isothermal. Inelasticity does not occur in these components as a result 
of mechanical loading alone but is induced mainly through thermal transient 

cycles and thermal gradients. Nevertheless, the constitutive equations 
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and damage models used in structural analysis and design are almost always 
based completely on experimental data collected under isothermal conditions. 
Isothermal tests are commonly conducted over the temperature range of 
interest; material "constants" are determined at each test temperature and 
then "fit" as functions of temperature across the temperature range to 
furnish a "noni sothermal " representation. This approach fails to reflect 
the strong thermomechanical path dependence observed in, for example, the 
cyclic hardening behavior of some alloys of interest - particularly in the 
presence of metallurgical changes (refs. 1-4). 

Not only are constitutive relationships used in large scale structur- 
al analysis often based solely on isothermal data but, further, they often 
represent only saturated cyclic behavior. As many common alloys continue to 
cyclically harden over several hundred or even several thousand cycles 
(ref. 1-4), and, as many large scale inelastic analyses are carried out only 
over a few cycles (largely for economic reasons), it is essential to charac- 
terize cyclic hardening over a wide range of cycling - not just for satura- 


ted conditions. 
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ABSTRACT 

Experimental evidence of' thermomechanical history dependence in the 
cyclic hardening behavior of some common high-temperature structural 
alloys is presented with special attention paid to the contribution of 
dynamic metallurgical changes. A discussion is given concerning the 
inadequacy of formulating "nonisothermal " constitutive equations sole- 
ly on the basis of isothermal testing. 

A representation of thermoviscoplasticity is proposed that qualita- 
tively accounts for the observed hereditary behavior. This is achieved 
by formulating the scalar evolutionary equation in an established vis- 
coplastic theory to reflect thermomechanical path dependence. Although 
the necessary nonisothermal tests for further quantifying the thermovis- 
coplastic model have been identified, such data are not yet available. 

To assess the importance of accounting for thermomechanical history 
dependence in practical structural analyses, two qualitative models are 
specified; the first is formulated as if based entirely on isothermal 
information, as is most often done; the second is made to reflect 
thermomechanical path dependence using the proposed thermoviscoplastic 
representation. Comparisons of predictions of the two models are made 
and a discussion is given of the impact the calculated differences in 
deformation behavior may have on subsequent lifetime predictions. 
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In this paper, we address thermomechanical path dependence in cyclic 
hardening with special emphasis on that induced through metallurgical chan- 
ges e.g., through dynamic strain aging. Definitive phenomenological evi- 
dence of thermomechanical history dependence is cited and a heuristic dis- 
cussion is given of its microscopic origin in solid solution hardening 
alloys. Related thermomechanical hereditary behavior is found in two-phase 
superalloys (ref. 5) arising from precipitation of the y' phase. 

A mathematical representation of thermoviscoplasticity is proposed that 
qualitatively accounts for the observed hereditary behavior. This is accomp- 
lished by appropriately formulating the evolutionary equation for the scalar 
state variable in a viscoplastic model developed principally by Robinson 
(refs. 6,7) to reflect thermomechanical path dependence. Although candi- 
date noni sothermal tests providing the necessary information to further 
quantify the thermoviscoplastic model have been identified, they will not be 
conducted until the presently expanding Structures Division Laboratory at 
NASA/Lewis Research Center becomes operative. 

To assess the importance of accounting for thermomechanical history in 
practical applications, two (qualitative) models are specified; one is taken 
as if based entirely on isothermal information as is traditionally done, and 
the second is taken to reflect thermomechanical path dependence using the 
representation proposed in the present work. Applications of each model are 
made to homogeneously stressed elements under some prescribed thermomechani- 
cal histories and the results are compared. 

THERMOMECHANICAL HISTORY DEPENDENCE AND METALLURGICAL CHANGE 
It is well known that cyclic hardening of some structural alloys 
within their temperature range of interest is influenced by the phenomenon 
of dynamic strain aging (ref 8). Strain aging occurs in solid solutions 
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where solute atoms are particularly free to diffuse through the parent 
lattice. It is energetically preferable for these solute atoms to occupy 
siles in the neighborhood of mobile dislocations where their presence 
immobilizes the dislocations or at least makes their movement difficult, 
thus causing strengthening. 

Isothermal cycling at temperatures where such metallurgical changes 
occur might therefore be expected to show abnormal hardening, i.e., higher 
hardening rates and greater saturation strengths than at temperatures both 
lower and higher. Evidence of strain aging in type 304 stainless steel is 
indicated in figure 1 (solid curves). There, the hardening rate and "sat- 
uration" strength are seen to be greatest at the intermediate temperature 
593C (1100F). Figure 2 shows the same behavior in the nickel-based alloy 
Hastelloy X (solid curves) at the specified strain range and rate and in 
the temperature range 427-649C (800-1200F). In this case the hardening 
rate and strength are greatest for the highest temperature 649C (1200F). 
For temperatures greater than 649C (not shown) both the hardening rate and 
saturation strength decrease markedly. Similar behavior is observed in 
type 316 stainless steel in figure 3. Further evidence of dynamic strain 
aging is reported in ref. 1-4. 

The observed hardening peak in these examples is interpreted as a 
manifestation of dynamic strain aging. At lower temperatures the mobility 
of solute atoms is far less and strain aging does not occur; at higher 
temperatures normal recovery processes, e.g., climb of edge dislocations, 
take over. 

In the aging process described dislocations can break away from their 
"Cottrell solute atmospheres" becoming mobile again. Although temporarily 
freed, dislocations can again be immobilized as solute atoms gradually 
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diffuse back to them. As the thermally activated process of diffusion is 
involved and solute atoms are migrating to dislocations which themselves 
are moving under the applied stress, it is not surprising that the ensuing 
inelastic deformation (cyclic hardening in particular) has a complex de- 
pendence on thermomechanical history. 

Phenomenological evidence of thermomechanical path dependence under 
cyclic conditions is seen in the results of the simple stepwise noniso- 
thermal tests reported in figure 1 and 2 (dotted curves). In these tests 
cycling is initiated at one temperature, after which the temperature is 
changed and cycling resumed. 

Figure 1 shows the result of a noni sothermal test on type 304 stain- 
less steel in which cycling is initiated at 427C (800F) for about fifty 
cycles; the temperature is then changed to 593C (1100F) and cycling con- 
tinued (dotted curve). Note that the material softens with the temperature 
increase, contrary to the implication of the isothermal hardening data. 
Hardening continues to 1000 cycles, however the stress amplitude remains 
measurably less than that indicated by the 593C isothermal curve. 

Figure 2 includes the results of four stepwise noni sothermal tests on 
Hastelloy X (dotted curves). Test path oabc represents the response to 
cycling at 649C (1200F) to 800 cycles and then at 427C (800F) to 2000 cy- 
cles. Test odef represents cycling at 649C (1200F) to 1600 cycles followed 
by cycling at 427C (800F) to 2000 cycles. Path oghj shows cycling to 1600 
cycles at 427C (800F) where the temperature is brought to 649C (1200F) and 
cycling continued. Finally, the two step test okmnps is reported. Here, 
cycling starts at 427C (800F) to 800 cycles, is changed to 649C (1200F) 
through 1200 cycles and then back to 427C (800F) to 2000 cycles. Note that 
the points f, c, s and 4 in figure 2 all represent the cyclic strength 
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(stress range) at 427C (800F) and 2000 cycles, the distinction being that 
each is reached over a different thermomechanical path. The range of 
"saturation" strengths at 4270 (800F) is close to 200 MPa. Similarly, 
points 1 and j in figure 2 correspond to 649C (1200F) strengths at 2000 
cycles following different thermomechanical histories. 

The features of these test results that reflect thermomechanical 
history dependence are: 

1) The change in strength (stress range) with temperature at a fixed 
number of cycles is always negative, i.e., an increase in temper- 
ature always results in a decrease in strength and vice versa - 
contrary to the implications of the isothermal data; 

2) The strength, in particular the "saturation" strength, depends on 
the temperature-strain history. The data suggest that the ma- 
terial retains a full memory of the thermomechanical path to 
cyclic saturation. Evidently, the information contained in the 
isothermal data is not sufficient to predict the hardening behav- 
ior under nonisothermal conditions. 
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REPRESENTATION OF THERMOVISCOPLASTIC CYCLIC HARDENING 


It is assumed in the following that the gradual hardening (or soften- 
ing) that accumulates over several cycles under cyclic stressing occurs 
isotropically and therefore can be described in terms of a single scalar 
state variable. That is not to say that kinematic (induced anisotropic) 
effects are not strongly present during cycling as the internal dislo- 
cation structure alters with immobilization and remobilization of dislo- 
cations under repeated stress reversals. Kinematic features are modeled 
through a tensorial state variable commonly identified as the internal or 
back stress. 

The flow law in the viscoplastic theory of refs. 6 and 7 is written 
(see Appendix) : 


■ f < F > 



a) 


in which e-jj denotes the components of inelastic strain-rate, p is a 
material parameter, f and F are material functions and J 2 is the second 
principal invariant of the effective stress 




( 2 ) 


as defined in the Appendix. 

For an initially isotropic material, the stress dependence enters 
through 


F = 



1 


(3) 
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F plays the role of a Bi ngham-Prager yield function in which k denotes 
the Bingham threshold shear stress below which the inelastic strain rate 
vanishes. In this stale variable description k is a scalar state variable 
accounting for isotropic hardening effects and a.jj are the components of a 
tensorial state variable (internal stress) accounting for kinematic 
effects. 

As the experimental work cited here is uniaxial we write equation (3) 
in uniaxial terms using 


thus 



(g-g ) 2 

3 

/ n2 

(q-g) 

K 


( 4 ) 

(5) 


where K=3k2. a and a denote the uniaxial components of the applied and 
internal (back) stress. Our concern in this study reduces to the speci- 
fication of an appropriate form of evolutionary law for K in equation (5) 
that reflects the thermomechanically hereditary behavior observed in figs. 
1-3 and cited in refs. 1-4. 

In the subsequent development we take isotropic hardening, through K, 
to depend on inelastic work 


W / Si jdei j , (6) 

however, the formulation is not innately limited to dependence on W and 
could, if appropriate, be taken to depend on some other scalar measure of 
mechanical deformation, say, 

P = | \jd e -j jde-j j (7) 

the accumulated plastic strain. 
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It is supposed that a typical material of concern cyclically hardens 
under some fixed strain range and strain rate as shown schematically in 
fig. 4. This is qualitatively in keeping with the actual material re- 
sponse observed in figs. 1-3 and with the physical arguments made in the 
previous section. Fig. 4a can be interpreted as a plot of Aa vs N (stress 
range vs number of cycles) or K vs W for isothermal cycling (solid curves) 
at the three temperatures Ti<T 2 <T 3 . Fig. 4b which is a plot of Aa or K vs 
temperature T shows the commonly observed "strain-aging peak". 

Classically, K has been taken as an explicit function of W and T, 
(e.g., ref. 9). This implies the existence of a unique surface K(W,T) as 
in fig. 4c that can be determined fully on the basis of isothermal harden- 
ing curves. This is equivalent to adopting the evolutionary equation for 
K as 

dK ' Iff d “ + It dT < 8 > 

which is a perfect differential and integrable independently of the ther- 
momechanical path W(T) as illustrated in fig. 4c. A candidate represen- 
tation of the isothermal curves of figs. 4a and 4b and the surface K(W,T) 
in fig. 4c is: 

-W/W (T) 

K(W,T) = K S (T) + [K-j (T) - K s (T)]e 0 (9) 

The functions K( 0 ,T) =K-j (T) and K(°°,T)=K S (T) are as illustrated in fig. 4b. 
The function W 0 (T) determines the temperature dependent cyclic hardening 
rate. 

The classical representation described above clearly does not reflect 
the thermomechanical path dependence apparent in the test data of figs. 

1-3 or in the schematically represented noni sothermal responses of fig. 4a 
(dotted curves). In particular, the saturated strength in the classical 
description depends only on the current temperature, i.e., K s (t ) 9 anc j j s 
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thus independent of the thermomechanical path leading to saturation. Actual 
data shows that the saturation strength does significantly depend on thermal - 
strain history. Furthermore, the classical description permits an increase 
in strength with an increase in temperature in the neighborhood ol the 
"strain-aging peak" (fig. 4b). Noni sothermal data shows that an increase in 
temperature is always accompanied by a decrease in strength, i.e., a decrease 
in K. 

A form of evolutionary equation for K in equation (5) that does allow 
for thermomechanical history dependence of the cyclic hardening response 
is the following: 

dK = r(W,T)dW + 0 (W,T)dT (10) 

Equation (10) represents a kinetic growth law for K which is not a perfect 
differential and is integrable only when the thermomechanical path W(T) is 
known*. The functions r(W,T) and 0(W,T) are independent and not related as 
3K/3VJ and 3K/3T are in equation (8). The function r(W,T) can be determined i 
the basis of isothermal testing whereas 0(W,T) cannot. A complete descrip- 
tion of the noni sothermal hardening behavior in terms of equation (10) 
necessarily involves nonisothermal testing. 

Candidate nonisothermal tests providing definitive information about th 
function 0(W,T) have been identified and are included in the test plan for 
the currently expanding Structures Division Laboratory at NASA Lewis Researc 
Center. For present purposes the results of related tests, identified ear- 
lier (refs. 1,2) by D. N. Robinson and conducted by R. W. Swindeman at Oak 
Ridge National Laboratory on types 304 and 316 stainless steel, have been us 

* It may be appropriate under some circumstances to account for thermal re- 
covery effects in the evolutionary law for K. As the emphasis here is on 
thermomechanical history effects, particularly in accompaniment of metallurg 
cal changes that occur at intermediate temperatures, we ignore thermal recov 
ery in K. 
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in suggesting an appropriate qualitative form of 0(W,T), and similarly of 
r(W,T) . On that basis we propose the following functional forms of r and 0 in 
the evolutionary equation (10): 

K (T) - K . (T) -W/W (T) 

r(W,T) = ( — - - g - p y 1 ) e ( 11 ) 


e(«,T> - - JWSL e- Q(H,(1/T » - 1/T) <1 2 > 

r 

where W is the inelastic work as defined in equation (6) and T is the 
absolute temperature. 

Note that r is identical to 3K/3W in equation (8) with K taken as in 
equation (9). Thus the functions K-j , K s and W 0 play the same roles as des- 
cribed earlier in the classical representation. However, the Arrhenius form 
adopted for 0 is not equal to 3K/3T in equation (8). The "activation energy" 
Q is, as yet, an unspecified function of W and T 0 is a reference temperature. 

For isothermal cycling at a fixed temperature T, we have dT = 0 and 
equation (10) integrates to equation (9). Thus, the classical representation 
where K in equation (5) is taken as an explicit function of W and T, (i.e., 
equation (9)), and the present formulation based on the evolutionary law ex- 
pressed in equations (10) - (12) give identical results under isothermal con- 
ditions. 

With temperature changes in the absence of mechanical deformation, i.e., 
with dW = 0, equation (10) gives 


= e(W,T) 


Q(W) -Q(W)(1/T 0 - 1/T) 

n 6 


(13) 


which, for Q(W)>0, ensures that 


dK 

dT 


< 0 


(14) 
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for all W and T - in keeping with the experimental observations cited above 
and contrary to the predictions based on the classical formulation. 

A finite temperature change from the reference temperature T 0 to T in 
the undeformed (W=0) state causes a change in K(=K-j) found by integrating 
equation (13) 




Q -Qo(l/T 0 - 1/T) 
o e 


(15) 


in which K 0 = Kj(T 0 ) and Q 0 = Q(0). Equation (15) integrates to 


Kid) = K 0 - (1 - 


e -Q 0 U/ T o- 1/T) j 


(16) 


The Arrhenius temperature dependence for Ki (initial strength) is substantia 
ted in biaxial (tension-torsion) experiments conducted on type 316 stainless 
steel (ref. 10). 

A critical feature of the evolutionary law expressed in equations (10)- 
(12) is that it permits the cyclic hardening, in particular the saturated 
strength, to be dependent on thermomechanical history. This feature will be 
explored in more detail in the following sections. 


SPECIFICATION OF TWO CYCLIC HARDENING MODELS FOR COMPARISON 
In this section we formulate two models for qualitatively describing t'r 
cyclic hardening behavior discussed above. Each is based on the viscoplasti 
constitutive equations outlined in the Appendix. The first. Model A, adopts 
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the "classical" representation of cyclic hardening identified in the preced- 
ing section and is treated as if based entirely on isothermal testing. The 
second model. Model B, is based on the thermoviscoplastic representation out- 
lined in the last section and is taken as being based on noni sothermal exper- 
imentation. Although, as indicated earlier, definitive nonisthermal data are 
not presently available, some qualitative guidance has been provided in for- 
mulating Model B from preliminary tests conducted at Oak Ridge National 
Laboratory. 

Predictions will be made using each model and compared in an effort to 
assess the importance of accounting for thermomechanical path dependence in 
practical applications. 

The details of the formulation of each of the models will now be given. 

1. Model A 

Constitutive model A assumes that the state variable K in equa- 
tion (5) is an explicit function of W and T, namely, that given by 
equation (9). 

The function K-j(T) in equation (9) is taken as specified in 
equation (16) with 


K 0 = 2 

Q 0 = 4000 (17) 

and T 0 = 800 K 


The functions K S (T) and W 0 (T) are, respectively 
K S (T) = A + B(T-Tp) + C(T - T p )2 


( 18 ) 
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with A = 16 

B = -0.04 
C = -0.001 

and W 0 (T) = A + B(T-T p )2 

with "A = 0.5 

F = 5xl0- 5 

In each case T p = 700 K . 

Isothermal curves K vs W for model A based on equation (9) and equation; 
(17) -(21) are shown in fig. 5a (solid curves). These represent hardening al 
6Q0K, 70QK and 800K. Fig. 5b shows the curves K(0,T) = K-j(T) and K(«,T) = 
K S (T). The curves in figs. 5a and 5b agree with the schematic representatio! 
in figs. 4a and 4b, including the appearance of the "strain aging peak". Noi 
that the surface K(W,T) defined by the solid curves in figs. 5a and 5b is 
unique in the case of Model A, and K remains on this surface for all thermo- 
mechanical paths as schematically illustrated in fig. 4c. Thus, Model A doe: 
not reflect thermomechanical path dependence. 

2. Model B 

Model B incorporates the kinetic growth law given in equations 

(10) - (12). The functions K-j(T), K S (T) and W 0 (T) in r(W,T)[equati 

(11) ] are taken identical to those defined for model A. Although a 
exact form of the "activation energy" Q(W) in equation (12) is not 
presently known, data in ref. (1,2) on type 304 stainless steel sug- 
gest it may not be a strong function of W. With that as a guide we 
choose Q to be a constant until appropriate test data become avail- 
able, i.e. we take Q(W) = Q(0) = Q 0 . The value of Q 0 has previous! 
been specified in equation (17). 


(19) 

( 20 ) 
(21) 
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The isothermal K vs W curves for model B are, of course, iden- 
tical to those for model A and likewise appear as the solid curves 
in figs. 5a and 5b. Now, however, thermomechanical paths other 
than isothermal ones do not necessarily remain on the surface de- 
fined by the isothermal data. This is illustrated by the calcula- 
ted curves (dotted lines) in fig. 5a depicting noni sothermal cycling 
based on Model B. These are in qualitative agreement with the 
experimental responses observed in figs. 1 and 2. 

The calculations made using either model A or B require coupling 
with a representation of linear thermoelasticity. As all calcula- 
tions are uniaxial and assume initial isotropy, it suffices to spe- 
cify only the Young's modulus E and the linear coefficient of ex- 
pansion a t . Thus: 


E(T) 

= 4x104 - 21.6 T 

(22) 

a t 

= 14 x lO" 6 

(23) 


giving E in ksi and a^- in K"l. Here only E is taken to depend on 
temperature. 

The two models formulated above, each based on the viscoplastic 
model stated in the Appendix, are valid within the following 
range of absolute temperature. 

600K < T < 800K. 

This completes the specification of the two models A and B. 

PREDICTIONS BASED ON MODELS A AND B 

We first consider predictions of uniaxial isothermal cycling over 
a fixed strain range and strain rate. Figures 6,7, and 8 show pre- 
dicted isothermal stress-strain loops for cycling over a strain 
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range of Ae = ±0.25% and a strain rate of e = 0.4 %/m at the temperatures 
600K, 700K and 8Q0K respectively. The hardening rate and saturation strength 
are evidently greatest for the intermediate temperature 70QK, as is typical 
of solid solution alloys that exhibit strain aging. The predicted isothermal 
hysteresis loops are identical for models A and B. A plot of stress ampli- 
tude Ao versus the number of stress reversals 2N corresponding to the iso- 
thermal hysteresis loops of figs. 6 - 8 is shown in figure 9 (solid curves). 

For thermomechanical paths other than isothermal ones, models A and B 
show qualitatively different behavior. For example, consider cycling at 
6Q0K over the path O^a in figure 9. At point a we suppose the temperature 
to be changed to 700K and cycling continued. The predicted response accor- 
ding to model A is an increase in A a to the point c at the next reversal and 
then to saturation at Sg under continued cycling. The stress amplitude S 2 
is precisely that at which saturation occurs for isothermal cycling at 700K. 
The prediction of model B, on the other hand, follows the path abSa. This 
shows a decrease in strength to b with the temperature increase, qualitative' 
ly in keeping with experimental observations, and saturation at S 4 , consider 
ably different than S 2 ; this too, is in accordance with experiments. Other 
predicted paths based on model B are 02 deSg and 03 ghSg. The former corres- 
ponds to cycling first at 700K ( 02 ^) then changing to 600K (eSs). The latte 
involves cycling at 800 K( 03 g) followed by cycling at 700K (hSg) with satura- 
tion occuring at Sg. Note that saturation at the temperature 700K occurs in 
these predictions based on model B at S 2 , S 4 or Sg depending upon the thermo 
mechanical path. 

We next consider thermomechanical cycling as illustrated in figure 10. 
Here, a uniaxial specimen is supposed to be initially at 6Q0K in its virgin 
(zero stress-zero strain) state. The total strain is held zero while 
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the temperature is cycled between 600K and 800K. Predictions of stress 
versus mechanical strain corresponding to a temperature ramp rate of 5K/sec 
are shown in figures 11, 12 and 13. Figure 11 shows the response based on 
inode! A; figure 12 shows that predicted from model B. Figure 13 is based 
on model A using only saturated isothermal data, i.e., taking K = K S (T) 
in Eq.(5), as is commonly done in practice. 

The limit cycles in figs. ll(or 13) and 12 are measurably different. 

The shape in figure 12 calculated using model B appears far more representa- 
tive of the results of generic experiments of this type (ref. 2). The 
difference in subsequent life predictions based on the two stable loops is 
speculative at present. Based on intuition developed for isothermal cycling 
one might expect the two lifetime predictions to be not too different since 
the inelastic strain range (width of the loops) and the stress values at the 
extremities of the loops are not drastically different. Nevertheless, the 
deformation prediction based on isothermal information alone (model A) does 
tend to indicate a narrower loop and thus might be expected to give the 
longer estimate of lifetime. Of course, such calculations depend strongly 
on the particular lifetime model being employed. 

Figures 14, 15 and 16 represent the same test as in the previous three 
figures except with the temperature ramp rate increased by an order of mag- 
nitude to 50K/sec. This calculation is made mainly to demonstrate the rate 
dependency of the viscoplastic constitutive model and the influence this has 
on thermomechanical response. Comparing the limit cycles in figure 14 or 16 
(model A) and figure 15 (model B) we see that their differences tend to be 
further exaggerated at the higher strain rate. Here, model A appears even 
more likely to give the less conservative estimate of lifetime. 
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The results of a final calculation are shown in figure 17. We again 
consider the problem illustrated in figure 10 but with the maximum tempera- 
lure taken as 700K instead of BOOK. The temperature rate is 5K/sec. In fig- 
ure 17 we compare the stable hysteresis loops based on model A, in which just 
saturated isothermal information is considered, and model B. The responses 
are now vastly different. Model A predicts shakedown (elastic cycling) 
about a sizable compressive mean stress. Model B does not show shakedown; 
its stable loop has finite width. Furthermore, the model B prediction indi- 
cates a far more damaging tensile mean stress. Here again, the trend is for 
the model based only on isothermal information (model A) to give the least 
nonconservative lifetime prediction. 

SUMMARY AND CONCLUSIONS 

Experimental evidence of thermomechanical history dependence in the 
cyclic hardening behavior of some structural alloys of interest is cited 
with particular emphasis on the contribution of dynamic metallurgical chan- 
ges. On this basis a discussion is given of the inadequacy of formulating 
"noni sothermal " constitutive relationships solely on the basis of isother- 
mal test data, as is often done in practice. 

A representation of thermoviscoplasticity is proposed that qualitative- 
ly accounts for the observed hereditary cyclic hardening behavior. This is 
accomplished by appropriately reformulating the evolutionary equation for 
the scalar state variable in an established viscoplastic theory to reflect 
thermomechanical path dependence. Although the necessary nonisothermal 
tests for further quantifying the thermoviscoplatic model have been tenta- 
tively identified, such data are not yet available. 
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To assess the importance of accounting for thermomechanical history 
dependence in practical structural analyses, two qualitative models are 
specified; the first is formulated as if based entirely on isothermal infor- 
mation, as is commonly the case; the second is made to reflect thermomechan- 
ical path dependence using the thermoviscoplastic representation proposed 
here. Comparisons of predictions of the two models are made and a specula- 
tive discussion is given of the impact the calculated differences in defor- 
mation behavior may have on subsequent lifetime predictions. 

Conclusions drawn from the present study are as follows: 

- Cyclic hardening occurs over a few thousand cycles in the alloys 
addressed here and for the strain ranges and rates considered. At 
lesser strain ranges hardening may continue over several thousand 
cycles. Therefore, a meaningful representation of cyclic defor- 
mation behavior should not be based solely on cyclically saturated 
strength, as is commonly the case. 

- The information contained in isothermal cyclic hardening data is not 
sufficient to predict hardening behavior under general thermomechan- 
ical conditions. 

- For the common alloys considered here, it appears that the satura- 
tion strength under cyclic conditions depends on the thermomechani- 
cal history to saturation, particularly in presence of metallurgi- 
cal changes. 

- The thermoviscoplastic representation presented here qualitatively 
accounts for the thermomechanical path dependence observed in the 
simple noni sothermal tests cited. 
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Comparison of predicted results of the two qualitative models indi- 
cated above under some simple thermomechanical loadings suggest that 
subsequent lifetime predictions may tend to be "nonconservative" 
when based on "noni sothermal " models formulated solely from isother- 


mal data. 
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APPENDIX 

Here we state the viscoplastic constitutive equations of ref. 6 and 7 
in noni sothermal and multiaxial form. As written, the equations are valid 
for "small" deformations and initially isotropic materials. 
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1 y y 

2 A i j A i j 


( 8A) 


I 0 = 4 a . .a . . 

2 2 ij ij 


( 9A) 


Here, e-jj denotes the components of inelastic strain rate, k and a^j are 
state variables and y, n, m, 3, R and H are material parameters. Speci- 
fication of the evolutionary law for the scalar state variable k is made in 
the foregoing text. k 0 is the value of k at the reference temperature T 0 . 

It is assumed that all of the temperature dependence in the flow law, equa- 
tions (1A), resides in the state variable k. Although this assumption is 
adequate for present purposes, a more generally valid representation might 
allow temperature dependence of the material parameters u and n, as well. 

The recovery parameter R is taken as being strongly temperature dependent as 
as indicated below. 

Values of the material parameters used in the present calculations are: 

y = 3.6 x 10 7 
n = 4.0 


3 = 0.75 
m = 3.87 

R = 8.97 x 10-8 exp[Q R (l/T 0 - 1/T) ] 

H = 9.92 x 10 3 

These are consistent with the units of ksi for stress, in/in for strain, 
time in hours and temperature in degrees Kelvin. The "activation energy" 

Q r = 4 x 10^ and the reference temperature T 0 = 800K. 

The values of the material parameters stated above were determined for 
a common steel (2 1/4 Cr-lMo), however, they are used here in the spirit of 
qualitatively representing a generic solid solution hardening alloy. 
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Figure 1 Isothermal cycling (solid lines) on type 304 stainless 
steel at 427, 593 and 700°C. Strain range is 0.54% and 
strain rate is 0.0001/s. Noni sothermal curve (dotted) 
shows cycling initially at 427°C changed to 593°C. 



stress range (MPa) 



Figure 2 Isothermal cycling (solid lines) on Hastelloy X at 427, 
538, 593 and 649°C. Strain range is 0.6% and strain 
rate is 0.001/s. Noni sothermal curves (dotted) obtained 
by cycling at 427°C and 649°C (see text). 
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Figure 3 Isothermal cycling on type 316 stainless steel at 350, 
450 and 550°C. Strain range is 0.56% and strain rate 





Figure 4 a) Typical isothermal (solid) and nonisothermal (dotted) 
hardening curves. 

b) Typical hardening vs. temperature curve showing strain 
aging peak. 

c) Surface k(W,T) 
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Figure 5 Isothermal hardening curves (solid lines) at 600, 700 
and 800 k for models A and B. Nonisothermal curves 
(dotted lines) corresponding to model B only. 
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Figure 8 


Isothermal hysteresis loop predicted by models A and B 
at 800 k. Strain range is 0.5%, strain rate is 0.004/m. 
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Fiqure 9 Isothermal cyclic hardening curves (solid) at 600 k . £9l5l ) 

700 k (0 o S 2 ) and 800 k (O 3 S 3 ) under conditions specified 

in figures 6 - 8 . Dotted curves depict nonisothermal responses 
of models A and B (see text). 
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Figure 12 Response predicted by model B for thermal cycling over 
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Figure 13 Response predicted by model A for thermal cycling over 
60:0 - 800 k at 5k/s including only saturated isothermal 
information . 
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Figure 16 Response predicted by model A for thermal cycling over 

600 -800 k at 50k/s including only isothermally saturate 
information . 
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